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Diffuse Gamma Flux from Antiproton Annihilation 

in our Galaxy 



The all-sky survey in high-energy gamma rays (E>30 MeV) carried out 
by the Energetic Gamma Ray Experiment Telescope (EGRET) aboard the 
Compton Gamma-Ray Observatory provides a unique opportunity to exam- 
ine in detail the diffuse gamma-ray emission. The observed diffuse emission 
has a Galactic component arising from cosmic-ray interactions with the local 
interstellar gas and radiation as well an almost uniformly distributed compo- 
nent that is generally believed to originate outside the Galaxy. The results 
of the observations have been interpreted as the extragalactic high-energy 
gamma-ray emission arising primarily from unresolved gamma-ray-emitting 
blazars. 

Here we consider another possible origin of the diffuse gamma-ray flux, 
namely, as originating from the annihilation of the antiprotons with the in- 
terstellar medium. 
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Abstract 



1 The annihilation cross sections 



The existing theoretical models based mainly on the partonic picture of the hadronic 
interactions are definitely invalid for pp annihilation at low energies. Due to this 
fact we used experimental data as for evaluation of the annihilation cross section as 
well as for the simulation of the final state and secondaries distribution. 

The annihilation cross section is the difference between total and inelastic ones 
c"ann ~ cTtot — o'ei- As it follows from data the dependence aann ~ is valid 
already for laboratory antiproton momenta piab < 1000 MeV/c. Thus, at Piab > 
300 MeV/c we used data from [|l| for the total and elastic cross sections and at 
momenta lesser than 300 MeV we used the dependence 



(7ann{P < 300 MeV) = do C{v*)/v* for annihilation cross section, 



The additional Coulomb factor C{v*) gives large increase for the annihilation 
cross section at small velocities of the antiproton and is defined by the expression 



where, Vc = ac, with a being the fine structure constant and c being the speed of 
light and v* is the velocity of the antiproton in the center-of-mass pp system. 



and 



a^i = const for elastic cross section 



i 




1 — exp (— 2 TT'Uc/'u*) 



We used the experimental data on the pp annihilation cross section 0, Q to 
normalize the 1/v behaviour and the Coulomb factor, choosing: 



^0 = <Z{P = 300 MeV) 



160 mb . 
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2 Gamma flux 



We used the experimental data |^ on the pp annihilation at rest to simulate the 
distribution of the final state particles, see Table. 



Table. Relative probabilities of pp annihilation channels. 



Channel 


Rel. prob., % 


Channel 


Rel. prob., % 




3 70 


2 7r+27r"r? 


60 




1.35 


/I p 


1.40 




1.35 




0.22 




9.30 


4.99 TT^ 


3.20 




23.30 




0.40 


7r+7r^47r'^ 


2.80 


2 7T+27T- 


6.90 




3.80 


3 7r^37r^ 


2.10 




7.30 


A'i^0.957r0 


6.82 


P'^TT^TT^TT^ 


3.20 




0.30 


P^TT^TT^TT^ 


3.20 




3.45 




16.60 




0.84 




4.20 




0.015 


3 tt+Stt^tt*^ 


1.30 




0.06 




1.20 







In practice, the approximation of the annihilation at rest is valid with very good 
accuracy up to laboratory momentum of the incoming antiprotons about 0.5 GeV 
because at this laboratory momentum the kinetic energy of the antiproton is still 
order of magnitude lesser than the twice antiproton mass. The average number of 
7's per annihilation is 

= 3.93 ± 0.24. 

The simulation of the distributions of final state particles has been performed 
according to phase space in the center-of-mass of the pp system. After that PYTHIA 
6.127 package p has been used to perform the decays of all unstable and quasi- 
stable particles. 

We used the spherical model for halo and parametrized the density distribution 
of interstellar hydrogen gas along the z direction (to North Pole) as: 



1 + {z/Df 

where n^^° = 10~^ cm~^ is the hydrogen number density in the halo, n^^'^ = 1 
cm~^ is the hydrogen number density in the disc and D = 100 pc is the half-width 
of the gaseous disc. The number of gammas arriving from the given direction is 
defined by the well known expression: 

L 



J,{E) = J dlij{E,r,z 



Here il^{E) is the density of gamma source along the observation direction / in 
suggestion of the isotropic distribution of the gamma source. The integration must 



be performed up to the edge of the halo L = —ax Rq + \J R\aio ~ -^o ~ ^x) with 
being the cosine of the line-of-sight to the x axis. 

The flux of the antiprotons has been chosen similar to the standard proton flux: 



Jp{Ekin) = Nfot /3 



, Ej^{ji J 



where (3 = p/E is the dimensionless antiproton velocity. The minimal kinetic energy 
of the antiprotons has been chosen equal to the velocity dispersion in the Galaxy, 
Vmin ~ 300 km/s, i.e.. 



EZ^ « 500 eV. 

Fig.|l](a,b) demonstrates the resulting differential gamma distribution in the Galac- 
tic North Pole direction in comparison with EGRET data |^ in the range 10 < 

< 1000 MeV. The peak of vr'' decay is clearly seen both in calculations as well as 
in experimental distributions. Fig.|^(c) shows the charged multiplicity distribution 
in the annihilation model described above. Experimental points have been taken 
from 

At this choice the necessary integral number density of the antiprotons Nfot is: 

Ntot = / 10"^^ cm-l 

Note here, that within reasonable variation of the index of the antiproton flux the 
integral number of the antiprotons does not varies significantly. This is due to the 
fact that the main contribution comes from the pp annihilation at the lowest velocity 
of the antiprotons. If we choose the minimal velocity of antiprotons according the 
velocity of the stellar wind, Vmnd ~ 1000 km/s, the necessary antiproton number 
density must be increased by order of magnitude and is Ntot ~ 10"^^ cm~^. This 
number also does not contradict to any theoretical or experimental limits on relative 
amount of the antimatter. 




Figure 1: Comparison of the calculated differential fluxes of 7 quanta from pp 
annihilation with experimental data EGRET on diffuse gamma background (a,b). 
There is also shown the charged multiplicity distribution in the annihilation model 
described in the text (c). 

3 Discussion. Antimatter Globular Cluster 

We considered the stationary case, provided by the presence of a permanent source 
of the antimatter. Assume [jlO| that antimatter stars in the Galaxy form globular 
cluster, moving along elliptical orbit in the halo. In the first approximation the 
integral effect we study depends on the total mass of the antimatter stars in the 
Galaxy and it does not depend on the amount of globular clusters. The main 
contribution into galactic gamma radiation comes from the antimatter lost by the 
cluster and spread over the Galaxy. 

There are two main mechanisms of antimatter loss by the cluster. The first one is 
the stationary mass loss by antimatter stars in the form of stellar wind. The second 
mechanism is the antimatter Supernova explosions. In both cases the antimatter is 
spread over the Galaxy in the form of positrons and antinuclei. Due to traveling 
in the halo magnetic fields the emitted antiprotons fill the whole halo and one can 
consider their number density to be constant over the halo. 

To evaluate the confinement time for the antiprotons in the Galaxy we used the 



results of the "two-zone" leaky box model [TT|. As the confinement time enters as 
a common factor in the predicted pp ratio [|11[|, we can find the necessary factor, 
performing the fit of the normalization to the experimental data. We removed from 
the fit two the most left points, strongly affected by the heliosphere [^]. Solid curve 
in Fig.0(a) presents the LBM predictions for the pp ratio, multiplied by the fitted 
factor K = 2.58, which factor increases the confinement time for slow antiprotons 
in the Galaxy up to 2 ■ 10^ years. Dashed curve is the phenomenological fit in the 
form R{E) = aE^~^^^^^. Experimental points have been taken from where 
references on the data can be found. 



R(E) = a*E**(b + c*lgE) 



a = 4.29E+00 ± 9.29E-01 
b = 1.52E+00 ± 3.1 1E-01 
c = -7.1 5E-01 ± 3.28E-01 
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Figure 2: (a) Fit of the pp ratio to experimental data. Solid line shows predictions 
of the two-zone leaky box model, increased by factor K 2.6. Dashed curve is the 
phenomenological fit, described in the text, (b) The comparison of the antiprotons 
annihilation time in the disc (solid) with the predicted LBM confinement time from 
Galaxy (dashed) and from disc (dot-dashed), (c) The free pass length of antiprotons 
in the disc. Dashed line shows the half-width of the disc. 



Fig.0 shows also the time of life of the antiprotons in the disc (b) and the free 
pass length of the antiprotons in the disc (c) versus antiproton kinetic energy. The 
correspondent curves for halo can be easily obtained if one reduces the hydrogen 
number density by two-three orders of magnitude according to H number den- 
sity in the halo. Note here that the presented in Fig.^ confinement time for very 



slow antiprotons has been extrapolated to small kinetic energies, suggesting that 
tesc{Ekin < 0.1 GeV^) ~ const = tesc{Ekin = 0.1 GcV) and is, in fact, the lower 
limit for the confinement time. 

From Fig.0(b) we can conclude that confinement time for low-energy antiprotons 
in the halo and in the gaseous disc is much greater than their life time with respect 
to the annihilation. Therefore we can consider the picture as a stationary one. From 
Fig.0(c) we see, that free pass length of the slowest antiprotons is lesser than the 
half-width of the disc and we can observe the products of the annihilation only but 
not the original antiprotons. The detailed analysis of the antimatter annihilation 
mechanism for Galactic diffuse gamma flux in the terms of antimatter globular 
cluster and the estimation of the total amount of antimatter stars and the expected 
antihelium flux for this case will be considered in the successive paper. 
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